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Abstract 

In this article, we consider the compressible Navier-Stokes equation with density depen- 
dent viscosity coefficients and a term of capillarity introduced by Coquel et al in |13| . 
This model includes at the same time the barotropic Navier-Stokes equations with vari- 
able viscosity coefficients, shallow- water system and the model of Rohde in [36] . 
We first study the well-posedness of the model in critical regularity spaces with respect 
to the scaling of the associated equations. In a functional setting as close as possible to 
the physical energy spaces, we prove global existence of solutions close to a stable equi- 
librium, and local in time existence for solutions with general initial data. Uniqueness is 
also obtained. 

1 Introduction 

This paper is devoted to the Cauchy problem for the compressible Navier-Stokes equation 
with viscosity coefficients depending on the density and with a capillary term coming from 
the works of Coquel, Rohde and theirs collaborators in [13], [36]. Let p and u denote the 
density and the velocity of a compressible viscous fluid. As usual, p is a non-negative 
function and u is a vector valued function defined on R N . Then, the Navier-Stokes 
equation for compressible fluids endowed with internal capillarity introduced in [36J reads: 

(SW) 



d t p + div(pu) = 0, 

d t (pu) + div(pu (g> u) - d\w(2p(p)Du) - V(A(p)divu) + VP(p) = npVD[p], 



supplemented by the initial condition: 

P/t=o = Po, pu/t=o = Po u o 

and: 

D[p] =<f>*p-p 

where <ft is chosen so that: 

cp € L°°(R N ) n C 1 ^) n W 1 ' 1 (R N ), [ 4>{x)dx = 1, (f) even, and cp > 
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and where P(p) denotes the pressure, p and A are the two Lame viscosity coefficients 
(they depend regularly on the density p) satisfying: 

H > 2p + NX > 0. 

(p is sometimes called the shear viscosity of the fluid, while A is usually referred to as 
the second viscosity coefficient). 

several physical models arise as a particular case of system (SW): 

• when k = (SW) represents compressible Navier-Stokes model with variable vis- 
cosity coefficients. 

• when k = and p(p) = p, X(p) = 0, P(p) = p 2 , N = 2 then (SW) describes the 
system of shallow- water. 

• when k = and p, A are constant, (SW) reduce to the Rohde model of chapter 
four. 

One of the major difficulty of compressible fluid mechanics is to deal with the vacuum. 
The problem of existence of global solution in time for Navier-Stokes equations was 
addressed in one dimension for smooth enough data by Kazhikov and Shelukin in |31j . 
and for discontinuous ones, but still with densities away from zero, by Serre in [38] and 
Hoff in [23]. Those results have been generalized to higher dimension by Matsumura 
and Nishida in [33] for smooth data close to equilibrium and by Hoff in the case of 
discontinuous data in [261 EZ] • 

Concerning large initial data, Lions showed in [32J the global existence of weak solutions 
for 7 > | for N = 2 and 7 > | for N = 3. Let us mention that Feireisl in [20] generalized 
the result to 7 > y in establishing that we can obtain renormalized solution without 
imposing that p G Lf oc , for this he introduces the concept of oscillation defect measure 
evaluating the loss of compactness. 

Other results provide the full range 7 > 1 under symmetries assumptions on the initial 
datum, see Jiang and Zhang [29|. All those results do not require to be far from the 
vacuum. However they rely strongly on the assumption that the viscosity coefficients are 
bounded below by a positive constant. This non physical assumption allows to get some 
estimates on the gradient of the velocity field. 

The main difficulty when dealing with vanishing viscosity coefficients on vacuum is that 
the velocity cannot even be defined when the density vanishes and so we cannot use some 
properties of parabolicity of the momentum equation, see [11], [12]. 
The first result handling this difficulty is due to Bresch, Desjardins and Lin in [7]. They 
show the existence of global weak solution for Korteweg system in choosing specific type 
of viscosity where p and A are linked. 

The result was later improved by Bresch and Desjardins in [4 J to include the case of 
vanishing capillarity (k = 0), but with an additional quadratic friction term rpu\u\ (see 
also [6]). However, those estimates are not enough to treat the case without capillarity 
and friction effects n = and r = (which corresponds to equation (1) with h(p) = p 
and g(p) = 0). 

The main difficulty, to prove the stability of (SW) , is to pass to the limit in the term 
pu®u (which requires the strong convergence of yj~pu) . Note that this is easy when the 
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viscosity coefficients are bounded below by a positive constant. On the other hand, the 
new bounds on the gradient of the density make the control of the pressure term far 
simpler than in the case of constant viscosity coefficients. 

In [6] Bresch and Desjardins show a result of global existence of weak solution for the 
non isothermal Navier-Stokes equation by imposing some condition between the viscosity 
coefficient and a bound by below on the viscosity coefficient. A. Mellet and A. Vasseur 
in using the same entropy inequality get a very interesting result of stability. They get 
more general estimate, which hold for any viscosity coefficients p(p), X(p) satisfying the 
relation: 

H(p) = p\'(p)-\(j>). (1.1) 

Mellet and Vasseur show in [33] the L 1 stability of weak solutions of Navier-Stokes com- 
pressible isotherm under some conditions on the viscosity coefficients (including (jl.ip ) 
but without any additional regularizing terms. The interest of this result is to consider 
conditions where the viscosity coefficients vanish on the vacuum set. It includes the case 
p(p) = p, X(p) = (when N = 2 and 7 = 2, where we recover the Saint- Venant model 
for Shallow water). The key to the proof is a new energy inequality on the velocity and 
a gain of integrability, which allows to pass to the limit. 

The existence and uniqueness of local classical solutions for (SW) with smooth initial data 
such that the density po is bounded and bounded away from zero (i.e., < p < po < M) 
has been stated by Nash in [35]. Let us emphasize that no stability condition was re- 
quired there. 

On the other hand, for small smooth perturbations of a stable equilibrium with constant 
positive density, global well-posedness has been proved in [33J. Many recent works have 
been devoted to the qualitative behavior of solutions for large time (see for example 
[244 I31j). Refined functional analysis has been used for the last decades, ranging from 
Sobolev, Besov, Lorentz and Triebel spaces to describe the regularity and long time be- 
havior of solutions to the compressible model [39\, [40] . [25] . [30] . The most important 
result on the system of Navier-Stokes compressible isothermal comes from R. Danchin in 
[15j and [18] who show the existence of global solution and uniqueness with initial data 
close from a equilibrium, and he has the same result in finite time. The interest is that 
he can work in critical Besov space (critical in the sense of the scaling of the equation.) 
We generalize the result of R. Danchin in considering general viscosity coefficient and 
in connecting this result with those of A. Mellet and A. Vasseur. This result improves 
too the case of strong solution for the shallow-water system, where W.Wang and C-J Xu 
in [H] have got global existence in time for small initial data with ho, uq £ H 2+s with 
s > 0. 

1.1 Notations and main results 

We will mainly consider the global well-posedness problem for initial data close enough 
to stable equilibria. Here we want to investigate the well-posedness of the system (SW) 
problem in critical spaces, that is, in spaces which are invariant by the scaling of the 
equations. Let us explain precisely the scaling of the system. We can easily check that, 
if (p,u) solves (SW), so does (p\,u\), where: 

p\(t,x) = p(X 2 t,Xx) and u\(t,x) = Xu(X 2 t, Xx) 
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provided the pressure law P has been changed into X 2 P. 

Definition 1.1 We say that a functional space is critical with respect to the scaling of 
the equation if the associated norm is invariant under the transformation: 

(p,u) — ► (p\,u\) 

(up to a constant independent of X). 

This suggests us to choose initial data (po,uq) in spaces whose norm is invariant for all 
A > by {p ,u ) — ► (po(A-), An (A-)). 

A natural candidate is the homogeneous Sobolev space 

fjN/2 x (i jAT/2-l ) JV ) but 

since 

fjN/2 j g no ^. i nc i uo ; ec [ i n l<x>^ we canno ^ expect to get L°° control on the density when 
P0 € H N ' 2 . 

This is the reason why as in the chapter two, instead of the classical homogeneous Sobolev 
space, we will consider homogeneous Besov spaces Bj? 2 x with the same 

derivative index. This allows to control the density from below and from above, without 
requiring more regularity on derivatives of p. In the sequel, we will need to control the 
vacuum, this motivates the following definition: 

Definition 1.2 Let p > 0, 9 > 0. We will note in the sequel: 



Let us first state a result of global existence and uniqueness of (SW) for initial data close 
to a equilibrium. 

Theorem 1.1 Let N > 2. Let p > be such that: P' (p) > 0, p(p) > and 2p(p) + 
\(p) > 0. There exist two positive constants So and M such that if qo G '~ , uq G 

B 2 x and: 

\\qo\\~N_ 1 N + IKH^w.! < e 

N 

then (SW) has a unique global solution (q, u) in which satisfies: 

\\(q,u)\\ E N < Af( lko|l 5 ^_i,^ + ||«o|l B ^_i), 
for some M independent of the initial data where: 

N , , , ~ N i N , i , j ~W i i AT,, r , iV ,, H 

= [C b (R + ,B~ 1 "T) nL 1 (M + ,S^ +1 'ir)] x [C b (R + , B^- Y ) N 

nL 1 (M + ,sf+ 1 ) iV ]. 

In the following theorems, we want to show some result of existence and uniqueness in 
finite time for large initial velocity and initial density close to some constant. 
The following result shows the existence and uniqueness in finite time with initial data 
in critical Besov space for the scaling of the equations. However as we said, we need for 
some technical reasons of an hypothesis of smallness on qo. We note that we work on 
Besov space Bp with general index p on the integrability. 
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Theorem 1.2 Let p E [1, +oo[. Let go £ -S/ an d u o S -B/ . Assume also that: 
\\Qo\\ < e /or a suitably small positive constant e > 0. 



R P 

p 



Under the assumptions of the theorem \l.l\ for the physical coefficients, there exists a time 
T > suc/i i/iai i/ie following results hold: 

N_ 

1. Existence: If p E [1,2N[ then system (SW) has a solution (q,u) in F p p with: 



F/ = C T (B P ) x (L\(B P ) n C T (B P 

JV 

2. Uniqueness: If in addition 1 < p < N then uniqueness holds in F p p . 
Moreover we have a control on the time T which may be bounded from below by: 

( , fJv_n IIA A-e- cVte2 \ eV 2 , 
mm 77, max (t > 0, > 2 Hy p ' A„u \\lp{ ~ )<~ — tt 

where v = min(/j(p), A(p) + 2/i(p)) and ?7o = |[«o]| js: • 

B p 

In the next theorem, we consider the case when the initial variational density go belongs 

■VP 



to p + B p p with e > and satisfies < p < p. Here we do not suppose a smallness 

condition on II go 1 1 & but we impose more regularity. 

B p 

Theorem 1.3 Let e E]0, 1[ and p E [1, and we assume that the physical coefficients 



-£ 1 

iV iV 



verify the same hypothesis as in theorem \l.l\ Let p$ E p + B p p for a constant p > 0, 

— +e— 1 

uq £ B p . Assume that there is a constant p such that: 

< p < po- 



There exists a time T > such that system (SW) has a unique solution (q, u) in F p + 
Moreover we have a control on the time T which may be bounded from below by: 



JV 



77,max(t > 0, V2 9( p" 1) ||A (? 77o||lp( = ) < ^r—fzr)). 

where: v = min(it(p), X(p) + 2p(p)) and Uo = \\uo\\ n . 

B p 

The present chapter is structured as follows. 

In the Section [2j we recall some basic facts about Littlewood-Paley decomposition and 
Besov spaces. 

In the Section [3] we prove the theorem 11.11 and in the Section 0] we show the theorem 11.21 
We conclude in the Section by the proof of the theorem 11.31 
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2 Littlewood-Paley theory and Besov spaces 
2.1 Littlewood-Paley decomposition 

Littlewood-Paley decomposition corresponds to a dyadic decomposition of the space in 
Fourier variables. 

We can use for instance any 92 G C°°(M N ), supported in C = {£ G R^/f < |£| < §} such 
that: 

Denoting /i = J 7-1 ^, we then define the dyadic blocks by: 

A ; u = (p(2' l D)u = 2 lN / h(2 l y)u(x - y)dy and Si« = V A fc u . 

Formally, one can write that: 

u = ^ A fc u. 



This decomposition is called homogeneous Littlewood-Paley decomposition. Let us ob- 
serve that the above formal equality does not hold in S (M. N ) for two reasons: 



1. The right hand-side does not necessarily converge in S 



2. Even if it does, the equality is not always true in S' (M N ) (consider the case of the 
polynomials). 

However, this equality holds true modulo polynomials hence homogeneous Besov spaces 
will be defined modulo the polynomials, according to [15] . 

2.2 Homogeneous Besov spaces and first properties 

Definition 2.3 For s G R, p G [l,+oo], q G [l,+oo], and u G S' (R N ) we set: 

ink, = (E( 2 'i a ^im 9 )*- 

lez 

A difficulty due to the choice of homogeneous spaces arises at this point. Indeed, ||.||b| 
cannot be a norm on {u G S (R N ), \\u\\b* < +00} because ||it||s» = means that u is 
a polynomial. This enforces us to adopt the following definition for homogeneous Besov 
spaces, see [15] . 



Definition 2.4 Let s G R, p G [I, +00], g G [1, +00]. 

Denote m = [s — ~] if s — ^ Z or g > 1 and m = s — — — 1 otherwise. 

• If m < 0, then we define Bp q as: 

B s p q = \u G S'{R N )/ \\u\\ B * tq < 00 and u = ^A ; u inS'^)! . 
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• If m > 0, we denote by V m \M> N ] the set of polynomials of degree less than or equal 
to m and we set: 

B s p q = !u £ S'(R N )/V m [^ N ] I IM|b», < oo and u = ^ A t u in S'(R N )V m [R N ] \ . 

The definition of B* r does not depend on the choice of the Littlewood-Paley decompo- 
sition. 

Remark 1 In the sequel, we will use only Besov space Bit with q = 1 and we will denote 
them by B* or even by B s if there is no ambiguity on the index p. 
Let us now state some basic properties for those Besov spaces. 

Proposition 2.1 The following properties holds: 

1. Density: If p < +oo and \s\ < N/p , then is dense in B^. 

2. Derivatives: there exists a constant universal C such that: 

C^IMIffi < [|Vu[| B «-i < C[|u[|s« . 

ii ii p,r ii n.Dp ir ii ii p,r 

3. Sobolev embeddings: If pi < P2 and r\ < r2 then Bp iri <— > B P2 ^ 2 n P2 . 

4- Algebraic properties: For s > 0, Bp r n -L°° is an algebra. Moreover, for any 

AT JV N_ 

p G [1, +oo] then B pl <— > B p p t00 n L°°, and B pl is an algebra if p is finite. 
5. Real interpolation: (B£ r ,B£ r ) e y = B^+(i-»)«a . 

2.3 Hybrid Besov spaces and Chemin-Lerner spaces 

Hybrid Besov spaces are functional spaces where regularity assumptions are different in 
low frequency and high frequency, see [To]. We are going to give the definition of this 
this new spaces and give some of their main properties. 

Definition 2.5 Let s, t G R. We set: 



\tl\ \ O s ;* 



■J2(2 qs \\\n\\ LP Yy + [^2(2^\\A q u\\ LP yy 

q<0 q>0 



Denote m = [s — ~] if s — ^ or r > 1 and m = s — ^ — 1 otherwise, we then define: 



• B s / = {u G S'(R N ) I \\u\lg.,* < +oo}, ifm < 

• B s / = {u£ S'{R N )/V m \R N ] I \u\^s,t < +00} ifm > 0. 



Let now give some properties of these hybrid spaces and some results on how they behave 
with respect to the product. The following results come directly from the paradifferential 
calculus. 
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Proposition 2.2 We give here some results of inclusion: 

1. We have B$ = B s pr . 

2. Ifs<t then Bpf r = B s pr n B pr or if s > t then B*;* = B s pr + B % pr . 

3. If si < s 2 and h > t 2 then Bp 1 / 1 ^ Bp 2 / 2 . 

Proposition 2.3 For all s, t > 0, 1 < r, p < +oo, i/ie following inequality holds true: 
||m;||s s ,t < C(||u||loo llull g s ,i + \\v \\u\\ „s,t ) . (2-2) 

For a// s\, s 2 , t\, t 2 < ^ such that min(si + s 2 , t\ + 1 2 ) > we have: 



TV 

ll'ut'llB 3 < C|M|_rs |H| jv if s < — . (2-4) 

For a proof of this proposition see |15j . The limit case s\ + s 2 = t\ + t 2 = in (12, 3D is 
of interest. When p > 2, the following estimate holds true whenever s is in the range 
(seee.fr [37]): 

lluull v < CIMIrs . (2.5) 

The study of non stationary PDE's requires space of type L p (0, T, X) for appropriate Ba- 
nach spaces X. In our case, we expect X to be a Besov space, so that it is natural to local- 
ize the equation through Littlewood-Payley decomposition. But, in doing so, we obtain 
bounds in spaces which are not type L p (0, T, X) (except if r = p) . We are now going to de- 
fine the spaces of Chemin-Lerner in which we will work (see [8] ) , which are a refinement of 
the spaces L^(B p r ). 

Definition 2.6 Let p G [l,+oo], T G [1,+co] and si,s 2 G M. We then denote: 

r T 



^2^(IIA/n(t)ir L p (iP) )" + C£ 2lrS2 ( f \\^iu(t)\\ p LP dty)y 



K0 l>0 

We note that thanks to Minkowsky inequality we have: 

IMIrP f5 s i' s 2'i ^ \\ u \\ tp ("R 3 i' 3 2\ if p^r. 
We then define the space: 

L p T (B s p ^ 2 ) = {u G L p t {B^ S2 )/\\u\\~ lPt {S;i ,s 2) < oo} . 

We denote moreover by Ct{Bp 1,S2 ) the set of those functions of L^?(i?p 1,S2 ) which are 
continuous from [0, T] to Bp 1 ' 32 . In the sequel we are going to give some properties of 
this spaces concerning the interpolation and their relationship with the heat equation. 
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Proposition 2.4 Let s, t, s±, s 2 E R, r, p, pi, p2 E [l,+oo]. We have: 
1. Interpolation: 



with l p = y 1 + k ^-> s = 9st + (1 - 0)s 2 and t = 0tx + [l- 6)t 2 . 
2. Embedding: 

L P T (B S /) ^ L P T (C ) and C T {sf) ^ C([0,T] x R d ) 

Here we recall a result of interpolation which explains the link of the space Bt, i with the 
homogeneous spaces , see [Tl"] . 

Proposition 2.5 There exists a constant C such that for all s E M, e > and 1 < p < 

+oo, we have 

» 1 j ST , , n ( ^ , 1 1 ^"1 1 _Sp ; 



\u\\b s <C rs 1 + log 



w Lbs 



To finish with we adapt the results of the paradifferential calculus on the product of 
Besov function to the spaces of Chemin-Lerner. So we have the following properties: 

Proposition 2.6 Letp,r E [l,+oo]. We have the two following properties: 

• Let s > 0, t > 0, l/p 2 + l/p 3 = 1/pi + 1/pa = l/p<l,ue L p T 3 (B S pfr) n L^(L°°) 
and v G L^{Bp*r) H Z^?(L°°). Then uv G ^(Bp.r) and we have: 

• Ifs x ,s 2 ,t 1 M < f , si + s 2 > 0, ii + i 2 > 1/ Pl + 1/ P2 = l/p<l,ue L^iB^ 1 ) 
dve L^(B s p f 2 ) then uv E L P T (B S 2 1+S2 ~ d/2 ) and 



an 



1 " vii L p t(b ;\ +s i-t .*i+*2-f } ~ II^IIl^^/ 1 ) 1 ^ 11 ^ 2 ^ 2 / 2 ) ' 



The analogous of the endpoint estimate (|2.5[) reads (for p > 2): 

_at < ||it||rPi/ RS JMIt^r-s v (2-6) 

whenever s is in the range (——, NN] and + = - < 1 (see the proof in |19|). For 
a proof of this proposition see |15j . Finally we need an estimate on the composition of 
functions in the spaces Lj,(Bp). 
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Proposition 2.7 Let s > 0, r G [l,+oo] and F G W^p 00 ^) suc/i i/tai F(0) = 0. 
There exists a function C depending only on s, p, N and F, and such that: 

W F ( U )h p T (B s p y 2 ) ^ C (\\ u \\l^(L^))\W\\ip t{ b s p ^)- 

Ifv,u£ L^iB'p 1 ' 82 ) n Lf(L°°) and G G W^ 3 ' 00 ^) then G{u) - G(v) belongs to 
Ly(Sp 1,S2 ) and it exists a constant C depending only of s,p,N and G such that: 

\\G{u) - G(v)\\ L p t {B s } ,s 2) < C(||u|Ug?(£«), \\v\\ L¥ - ^IIl^KV" 2 ) 

(! + IMk°?(X°°) + IMU§?(X°°)) + lb - u \\l™ (l°°)(\\ u \\lp, (Bpi; 32 ) + IMIl^V 52 ))- 

The proof is a adaptation of a theorem by J.Y. Chemin and H. Bahouri in pp. 

We end this section by recalling some estimates in Besov spaces for transport and heat 

equations. For more details, the reader is referred to [8] and [17j . 

Proposition 2.8 Let (p,r) G [1, +oo] 2 and s G (— min(~ 3 4-), y + 1)- Let u be a vec- 

JV 

tor field such that Vu belongs to L l (Q,T;Bp )r n L°°). Suppose that qo G Bp T , F G 
L 1 (0, T, Bp r ) and that q G L?p{Bp r ) D C([0,T];<S ) solves the following transport equa- 
tion: 

(d t q + u-Vq = F, 
\ q t =o = qo- 

Let U(t) = $1 ||Vn(r)|| n dr. There exits a constant C depending only on s, p and 

B» r C\L°° 

N , and such that for all t G [0, T\, the following inequality holds: 

kh rW < exp^W^golb^ + [^~ CU{T) \\nr)\\B lr dr^j 
If r < +oo then q belongs to C([0, T]; Bp r ). 

Actually, in [T7], the proposition below is proved for non-homogeneous Besov spaces. 
The adaptation to homogeneous spaces is straightforward. Let us now some estimates 
for the heat equation: 

Proposition 2.9 Let s G R, (p,r) G [l,+oo] 2 and 1 < p% < p\ < +oo. Assume that 
uq G Bp r and f G L^(Bp :r 2+2 ^ P2 ). Let u be a solution of: 

j d t u - ij.Au = f 
[ u t =o = uo ■ 

Then there exists C > depending only on N, p, p\ and p% such that: 

ihi^^M) ^ c (\\ u °H, r +/^~ 1 imizp 2(5r ^/p 2) ) • 

If in addition r is finite then u belongs to C([0, T], Bp r ). 
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The proof of local well-posedness for initial density bounded away from zero requires 
estimates in Bp spaces for the following linear system: 

{dtu — /2div(aVu) — (A + /Z)V(adivu) = G, 
u t =o = uo, 

where u(t, x) 6 ]R , v = 2/2 + A > and the diffusion coefficient is assumed to satisfy: 

< a < a(t,x) < a. (2.7) 

We can prove that the solution of the previous system satisfy estimates analogous to 
those of proposition 12.91 see pj)] . 

Proposition 2.10 Let 1 < p < +oo, 1 < ct\ < r < +oo and s be such that max(l, ~) < 
s < — + 1. Set a 2 = 2 n an d let ai be such that = + \. Let assumption \2. 7| j 

s v 012 

be fulfilled and u be a solution of \2. 3\) . We suppose that the regularity index r satisfies: 

2 N 2 

1 <r<l hs. 

a>i p a>i 

Then the following estimate holds for all a £ [r, +oo] : 
1 1 ix- 1 1 T+ 2 < \\uo\\b t + \\G\\~ r+z-2 + ll^ a ll ' 

3 Proof of theorem 11.11 
3.1 Sketch of the Proof 

In this section, we give the sketch of the proof of theorem 11.11 on the global existence 
result with small initial data. 

We will suppose that p is close to a constant state p, so that p will be strictly superior to 
a positive constant, we will use the parabolicity of the momentum equation to get a gain 
of derivatives on the velocity u. The density p has a behavior similar to the solution of a 
transport equation. Let us rewrite the (SW) system in a non conservative form in using 
the definition 11.21 



(SW1) < 



dtq + uSJq + divu = F 

p p 



— np4> * Vg = G 
where we have: 

F = —qdivu, 

G = A(p,u)+K(p)Vq, 



11 



where we set: 



A{ P ,u) = [ 



div(/x(p)£>(«)) 



K{p) 



pP'(p) 
P 



P 

-P'(P)- 



A»(p) A i , r V((M(p)+A(p))div«) n{p)+\{p) 
-^—Au\ + [ = 



Vdivu] , 



For s £ M, we denote A s /i = J r_1 (|^| s /i). We set now: 

d = A^divu and f2 = A -1 cur In 

where d represents the compressible part of the velocity and O the incompressible part. 
We rewrite now the system (SW1) in using these previous notations on a linear form: 



(SW2) 



where we have 



'd t q + Ad = F 1 , 
d t d - uAd - 5Aq + RA(4> * q) 

o t n - p,An = H x 

u - 



G^ 



A -1 Vd - Adivft 



P(p) 



P 



p 



2/i + X, 5 = K,p + P(p), and K = np. 



We have in our case: 



F\ = —q divu — u • Vq, 
d = -A- 1 div(G), 
Hx = -A _1 curl(G). 

The first idea will be to study the linear system associated to (SW2). We concentrate 
on the first two equations because the third equation is just a heat equation with a non 
linear term. The system we want to study reads: 

[d t q + Ad = F, 

\ d t d - vAd - 5Aq + RA(cp * q) = G' . 

This system has been studied by D. Hoff and K. Zumbrum in [28] in the case R = 0. 
There, they investigate the decay estimates, and exhibit the parabolic smoothing effect 
on d and on the low frequencies of q, and a damping effect on the high frequencies of q. 
The problem is that if we focus on this linear system, it appears impossible to control 
the term of convection u ■ Vg which is one derivative less regular than q. Hence we shall 
include the convection term in the linear system. We thus have to study: 

J d t q + v ■ Vq + Ad = F, 
\ d t d + v-Vd- vAd - 5Aq + RA(cf) * q) = G, 



(SW2) 



where v is a function and we will precise its regularity in the next proposition. System 

(SW2) has been studied in the case where <j> = by R. Danchin in [15j . we then adapt 

the proof in taking into consideration the term coming from the capillarity. 

In the sequel we will assume v > and 5 — RW^Wl 00 > c > 0. We get then the following 

proposition. 
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Proposition 3.11 Let (q,d) a solution of the system (SW2)' on [0, T[ , 1 — y < s < 



1 + y an d = Iq ll u ( T )ll n +1 dr. We have then the following estimate: 



ll(?»«Ollfl.-l,. xB -l+ / H(9> rf )( T )ll5H-l,. x fl.+l*" 

•/o 

<C7e cv W(||((jb > db)|| 5 ._ 1 ,. xfl ^ 1 + / t e- cv M||(J?,G)(r)|| 5 ._ 1 ,. xfl ._ 1 cfa-), 

JO 

where C depends only on v, 5, k, <j), N and s. 
Proof of Proposition 13. lit 

Let (q, u) be a solution of (SW2) and we set: 

q = e- KV ®q, u = e~ KV W U , F = e ~ KV ^F and G = e~ KV ®G. (3.8) 

We are going to separate the case of the low and high frequencies, which have a different 
behavior concerning the control of the derivative index for the Besov spaces. 
In this goal we will consider the two different expressions in low and high frequencies 
where Iq G Z, A, B and K\ will be fixed later in the proof: 



ff = Hmh-~<mA*qi) + \\di\\l2-2K 1 (Aq u d l ) for I < l Q , 
fl = \\Aqi\\h+Mdl\\% ~ ~(A®,d,) for I > Iq. 



(3.9) 



In the first two steps, we show that K\ and A may be chosen such that: 

tf(.-i) f 2 ^ 2 ^max(l,2- i )||^|l! 2 + &- X \&ifv> (3-10) 
and we will show the following inequality: 



\~<E^ amm (^ 2 '' ^f'l — C2~ l ( s ~ 1 ^a>ifi( \\(F, G)\\ £, s -i.s xB s-i 



+ V'\\(q,d)\\ Sa . 1 ^ Ba ^)-KV'ff. 



(3.11) 



where Yli^z a i — 1 ana ^ o is a positive constant. 

This inequality enables us to get a decay for q and d which will be used to show a 
smoothing parabolic effect on d. 

Case of low frequencies 

Applying operator A; to the system (SW2)' , we obtain then in setting: 

qi = Aiq, di = Aid. 

the following system: 

r ^-qi + A[(v • V§) + Adi = F- KV\t)q u 
d (3.12) 

L — + A t (v ■ Vdi) - vAd x - 5Aqi + kA(0 * q t ) =G X - KV\t)d x . 
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We set: 

fl=SMh + \\<k\\ 2 v >-2K 1 (A%,d l ) (3.13) 

for some K\ > to be fixed hereafter and (•, •) noting the L 2 inner product. 
To begin with, we consider the case where F = G = 0, v = and K = 0. Taking the L 2 
scalar product of the first equation of (|3. 12|) with q~\ and of the second equation with di, 
we get the following two identities: 



^^llftllia + (Adi.gi) = 0, 

T^INlis + *||M|lia " *(A<B,dl) + K(A(0 * qi), di) = 0. 



(3.14) 



In the same way we have: 

^(«,«*^) + (M,0*«) = O, (3.15) 
because we have by the theorem of Plancherel: 

We want now get an equality involving ^(Aci;, qi). To achieve it, we apply z^A to the first 
equation of (13.120 and take the L 2 -scalar product with di, then take the scalar product 
of the second equation with Aqi and sum both equalities, which yields: 

^(A^dO + HMIlia -m qi \\ 2 L2 +R\\<t>*Aq l \\ 2 L2 + P(A 2 d h A gi ) = 0. (3.16) 



By linear combination of (|3.14j) and (|3.16j) . we get: 

\j t ff + {v- JrOHMIIla + ^i(^||A^||| 2 - * A qi \\\ 2 ) - uK 1 {A 2 d h A qi ) = 0. 

(3.17) 

And as we have assumed that: 5 — k||</>||l°° > c > we get: 

\j t ff + (u - ^i)||A^||| 2 + K lC \\ qi \\l 2 - vKi(A 2 di,Aqi) < 0. (3.18) 

Using spectral localization for d\ and convex inequalities, we find for every a > 0: 

n2 2l ° 1 
|(A 2 d z ,A % )|<^-||A^||! 2 + -||A ft || 2 i2 . 

In using the previous inequality and (I3.17p . we get: 

--/ 2 + (i> — K\ - ^-)||A^|| 2 2 + (K lC - ^MQiWh < 0. (3.19) 
From (|3.13p and (|3. 19|) we get in choosing a = v and K\ < mm(-^-, -2 )1 then: 

+ a2 2l f 2 < 0, (3.20) 
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for a constant a depending only on v and K\. 

In the general case where F, G, K and v are not zero, we have: 

—ff + (a2 2l +KV')ff < {F h q{) + (G h d t ) - K{AF u di) - K(AG hqi ) - (A,(* • V§),^ 

- (A t (v • Vd),di) + K((AAi(v -Vq),di) + ((AA,(u • W), ®). 

Now we can use a lemma of harmonic analysis in [TB] to estimate the last terms, and get 
the existence of a sequence (a;)/ e z such that X^ez a ' — 1 an< ^- 

—ff + (a2* + KV')ff < atffi-^iW^G)^^^ + V'\\(q, ^lla-w^-i) . 

(3.21) 

Case of high frequencies 

We consider now the case where I > Iq + 1 and we recall that: 

f? = \\m\\h+MMh-=®A 

For the sake of simplicity, we suppose here that F = G = 0, v = and K = 0. We now 
want a control ||A<//||? 2 on e apply the operator A to the first equation of (|3,12|) . multiply 
by Aqi and integrate over M. N , so we obtain: 

^I|A<Z/Hi2 + (A 2 ^,A % )=0. (3.22) 

Moreover we have: 

IjMWh + - Si^di) + S(A(0 * q t ), d t ) = 0. 

(3.23) 

— (A^,^) + ||A^||| 2 -5||A^||| 2 +k||0*A%||2 2 + u(A 2 d h A qi ) =0. 
By linear combination of (|3.22|) - (|3.23|) we have: 



2^/l + ^II A «II2 2 + " r) [IMI& - ^(A®, d,) + AR(A{cp * ffI ). = 0- (3-24) 
Moreover we have: 

| - A5(A qil di) + AR{A((f> * q{),di)\ < A{5 + k||0|| l ~)|(A®, d t )\ 
We have now in using Young inequalities for all a > 0: 

|(^A ft )|<|||Aft||! 2 + ^||^|li 2 , 

So we get: 

\U + ^ A > ~ \ ~ i^' + { \ ~ 5 (3 ' 25) 
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So in choosing: 

1 2 

a = — 7 and A > max(— , 1) 
vA v 

there exists a constant a such that for I > Iq + 1 we have: 

In the general case where F, G, H, K and v are not necessarily zero, we use a lemma of 
harmonic analysis in [15] to control the convection terms. We finally get: 



\±ff + (a + KV')ff < aJfi-^iUF, G) ||g s _ lj Bs _ x 

i at ^3_ 27 ^ 

+ ^||(M|| g ._ 1 ,. xfl ._ 1 ). 

This finish the proof of $Bfy and <j!TTT]> . 
The damping effect 

We are now going to show that inequality (|3.11|) entails a decay for q and d. In fact we 
get a parabolic decay for d, while q has a behavior similar to a transport equation. 
Using hf = ff + 5 2 , integrating over [0, t] and then having 5 tend to 0, we infer: 

/,(t) + amin(2 2 ',l) /* /,(r)dr 

JO 

</ i (0) + C2-^- 1 ) /*aj(r)||(F(r) ) G(r))|| 5 ^ 1 ,. xfl .dr (3.28) 



o 



Thanks to (|3.1U|) . we have in taking i"T large enough : 

X)(C2- , t- 1 )a,(r)||(?,^|| 5 ._ 1 ,. xB .-Jir/ J (r)) <0, 



In multiplying (|3T28j) by 2 / ( s_1 ) and in using the last inequality, we conclude after sum- 
mation on Z, that: 

||g(i)b s -i, s + Wss-i+a /VWIIS.-M^ + E r«2^ s - 1 )min(2 2i ,l)||^(r)|U 2 dT 



/•t 



< H^,do)|lB.-i,»xB-i + / IIC F > G! )llB-i.«xB-i dr - 







(3.29) 
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The smoothing effect 



Once stated the damping effect for q, it is easy to get the smoothing effect on d by 
considering the last two equations where the term Aq is considered as a source term . 
Thanks to (|3.29j) . it suffices to prove it for high frequencies only. We therefore suppose 
in this subsection that I > Iq for a Iq big enough. 
We set gi = \\di\\i2 and in using the previous inequalities, we have: 

\j t \\di\\h + ^IIMIlls " S(Am, d t ) + R(A(<i> * qM) =G x -d x - KV'^UWli. 
We get finally with a > 0: 

~gf + a2 2l gf < ^(||Ag z || L 2 + ||G,|| £a ) + g l V\t)[Ca l 2~ l ^ 1 )\\d\\ BS -, - K 9l ). 

We therefore get in using standard computations: 

V^-^I^IUa+a [ t Y,2 l{S+1) \\dl(T)\\ L idT<\\d \\ B s- 1 + f \\G(r)\\ B s-rdT 
l>lo Jo l>l Jo 

+ fY,1 lS \\UT)\\v+CV{t) sup ([|d(r)|| fl .-x). 

JO ^ re[0,t] 

Using the above inequality and (|3.29p . we have: 

r^2^+ 1 )||^(r)|| L2( ir<(l + y(t))(||go||^- 1 , 3 + ||dob-) 
Jo l>lo 

(3.30) 



i7(r)|| g ._ 1 ,. + ||G(T)|| B .- 1 )dr. 



o 



Combining that last inequality (|3.30p with (|3.29p . we achieve the proof of proposition 

Em 

□ 

3.2 Proof of theorem 11.11 

This section is devoted to the proof of the theorem 11.11 The principle of the proof is a 
very classical one. We want to construct a sequence (q n ,u n ) n ^n of approximate solutions 
of the system (SW), and we will use the proposition 13.111 to get some uniform bounds 
on (g n ,'u n ) ne N- We will conclude by stating some properties of compactness, which will 
guarantee that up to an extraction, (q n , u n ) ne ^ converges to a solution (q, u) of the system 
(SW). 

First step: Building the sequence (q n ,u n ) n ^ 

We start with the construction of the sequence (q n ,-u n ) ng N 5 in this goal we use the 
Friedrichs operators (J n )neN defined by: 

Jn9 = ^ 1 (l B (± in) ?), 
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where T 1 is the inverse Fourier transform. Let us consider the approximate system: 



{ d t q n + J n (J n u n ■ VJ n q n ) + AJ n d n = F n 

d t d n + J n (J n u n ■ VJ n d n ) - vAJ n d n - 5AJ n q n -R<f>* AJ n q n = G n 

d t n n - HAJ n n n = H n 

u n = -A -1 V<T - A _1 divfi n 
k (q n ,d n ,Q n )/ t=0 = (J n qo, J n do, J n Qo) 



(3.31) 



with: 



F n = -J n ((J n q n )divJ n u n ), 

G n = J n A- 1 div[^l(^(p(l + J n q n )), J n u n ) + K(<p(p(l + J n q n ))Vq n ], 
H n = J n A- 1 curl[^l(^(p(l + J n q n )),J n u n ) + K{ V {p{\ + J n q n ))Vq n ] . 

where <p is a smooth function verifying (f(s) = s for - < s < n and tp >\. 
We want to show that (|3.3ip is only an ordinary differential equation in L 2 x 1? x L 2 . 
We can observe easily that all the source term in (|3.3ip turn out to be continuous in 
L 2 x L 2 x L 2 . As a example, we consider the term J n A(ip(p(l + J n q n )), J n u n ). We have 
then by Plancherel theorem: 



J n 



AW(p(tp(p(l + J n q n ))DJ n u n ) 
<p{p(l + J n q n )) 



)\\ L 2 < n\\fi(ip{p(l + J n q n ))DJ n u n \\ L 2 



l ip(p(l + J n qn)) 



\L° a ; 



< 4M n n 2 \\u n \\ L 2. 



where M n = \\p(ip(p(l + J n q n )\\L°°- 

According to the Cauchy-Lipschitz theorem, a unique maximal solution exists in C([0, T n ); L 2 ) 
with T n > 0. Moreover, since J n = J 2 we show that {J n q n , J n d n , J n Q n ) is also a solu- 
tion and then by uniqueness we get that (J n q n ,J n u n ) = (q n ,u n ). This implies that 
(q n , d n , J7 n ) is solution of the following system: 



and: 



' d t q n + J n {u n .Vq n ) + Ad n = Ff 
d t d n + J n (u n .Vd n ) - vAd n - 5Aq n -R(j)* Aq n = G\ 
d t Q n - uASl n = H? 
u n = -A^VcT - A _1 divfi n 

. W > d n , Q n )/t=o = ( J n qO: Jndo, Jn&o) 



F? = -J n (q n divu n ), 

G\ = J„A- 1 div[^l( 93 (p(l + q n )) , u n ) + K(<p(p(l + q n ))] , 
HI = J n A- 1 curl[^( V 9(p(l + q n )), u n ) + K(<p(p(l + q n ))] ■ 



(3.32) 
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And the system (|3.32p is again an ordinary differential equation in L 2 with: 

L 2 n = {g G L 2 (R Ar )/supp?C n)}. 

Due to the Cauchy-Lipschitz theorem again, a unique maximal solution exists in C 1 ([0, T^);L 
with T' n >T n > 0. 

Second step: Uniform estimates 

In this part, we want to get uniform estimates independent of T on ||(g n ,?/ n )|| jv for all 
T < T' n . This will show that T' n = +00 by Cauchy-Lipchitz because the norms || • || jv 

JV 

and I? are equivalent on L\. E^ Let us set: 
E(0) = \\q \\ s N_ h N + ||uo|| B £, 



E(q,u,t) = \\q\\ ~n , n + llgll jv , + llgll , ~jv jv + 



\ 



and: 

T n = sup{i G [0,T^),^(g n ,u n ,i) < 3C£(0)} 

C corresponds to the constant in the proposition 13.111 and as C > 1 we have 3C > 1 so 
by continuity we have T n > 0. 

We are going to prove that T n = T n for all n G N and we will conclude that Vra G N 
T n = +00. To achieve it, one can use the proposition 13, 1 II to the system (|3.32p to obtain 
uniform bounds, so we get in setting V n (t) = \\u n \\ x n +1 : 

||(<f> n )|| £f <Ce cv ^{\\q4^_ lif + \\uo\\ Bf + ^ V c ^)(||Ff (r)|| ~ f _ lif 

+ ||G?(r)|| B ^_ 1 + ||flr(r)|| fl4 _ l )dr.) 

Therefore, it is only a matter of proving appropriate estimates for i 7 ™, Gf and H™ in 
using properties of continuity on the paraproduct. 

We estimate now \\F™\\ 1 ~jv_ 1 jv in using proposition 12.61 and 12.71 

||F-, n || , jv , jv < C||<7 n l| jv , jv ||divu n || , jv , 

We now want to estimate G": 

\\A(tp(p(l+q n )),u n )\\ . jv x <C\\u n \\ . «+A\q n \\ , n ,(l + ||g n || , n), 

We can verify that K fulfills the hypothesis of the proposition 12.71 so we § e t ; 
\\K{<p{p(l + q n ))Vq n \\ . jv 1 <C||g n || 2 , ,v ||cf|| ~jv 
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Moreover we recall that according to proposition 12. 6t 

IU n ll 2 \\r, n \\ \\n n \\ 

" "l4(.BT) - UH "^(BT-W) 111 ll Ll(B1T+ 1 'V) 

We proceed similarly to estimate \\H™\\ n _ 1 and finally we have: 

l t (b^ ) 

PTII , , # 1 , + ||G?|| ., n u + ||flT|| iv ! <2C{E 2 {q n ,u n ,T) 

+ i5V,u n ,T)), 

whence: 

||(g n ,u n )|| w <Ce c23E(0) E(0)(l + 18CE{0)(l + 3E{0))), 
We want now to get: 

e 3C 2 £(0) (1 + 18C r S ( )(i + 3£(0))) < 2 

for this it suffices choose £?(0) small enough, let E(0) < £ such that: 

1 + 18CE(0)(1 + 3£(0)) < § and e 3c2£ (°) < 1 

So we get T n = T' n , indeed we have shown that VT such that T < T n : 

E(q n ,u n ,T) < 2CE(0). 

Then we have T n = T' n , because if T n < T' n we have seen that E(q n ,u n ,T n ) < 2CE(0) and 
so by continuity for T n +e with e small enough we obtain again E(q n ,u n ,T n +e) < 3CE(0) 
and stands in contradiction with the definition of T n . 
So if T n = T n < +00 we have seen that: 

E(q n ,u n X)<3CE(0). 

As ||o n || ~jv < +00 and ||u n || ~n , < +00, it implies that ||a n ||r°° it2\ < +00 

IIWII^oo II nn L™ (BT- 1 ) ' 1 » w » L T ^n) 

and ||it n ||^oc < +00, so by Cauchy-Lipschitz theorem, one may continue the solution 

beyond T n which contradicts the definition of T n . 

Finally the approximate solution (q n ,u n ) n( =N is global in time. 

Second step: existence of a solution 

In this part, we shall show that, up to an extraction, the sequence (q n ,u n ) n< =n converges 
in T> (R + x R ) to a solution (q, u) of (SW) which has the desired regularity properties. 
The proof lies on compactness arguments. To start with, we show that the time first 
derivative of (q n ,u n ) is uniformly bounded in appropriate spaces. This enables us to 
apply Ascoli's theorem and get the existence of a limit (q, u) for a subsequence. Now, 
the uniform bounds of the previous part provide us with additional regularity and con- 
vergence properties so that we may pass to the limit in the system. 
It is convenient to split (q n ,u n ) into the solution of a linear system with initial data 
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(q n ,u n ) and forcing term, and the discrepancy to that solution. 
More precisely, we denote by , u\ ) the solution to: 

d t ql + divu£ = 

d t ul - Au\ + Vql = (3.33) 

(QLi v l)/t=o = (J n qo,JnU ) 

where: A = p,A + (A + / u)Vdiv and we set (if 1 , u n ) = (q n — q^, u n — u^). 
Obviously, the definition of (q£, v^) /t=o entails: 

— JV JV 

(9l)/*=o ^ go in 5 2 A » 2 , (n2)/ t =o u in £ 2 \ 

The proposition 12.91 insures that (g£,«2) converges to the solution (qL,UL,) of the linear 
system associated to (|3.33p in E~ . We now have to prove the convergence of {q a ,u n ). 
This is of course a trifle more difficult and requires compactness results. Let us first state 
the following lemma. 

Lemma 1 (q n , u n )) ne n is uniformly bounded in (R + ; B^^ 1 ) x (C^(R + ;B^~^)) N . 
Proof: 



In all the proof, we will note u.b for uniformly bounded. 

We first prove that ^q n is u.b in L 2 (R + , B~~ 1 ), which yields the desired result for q n . 
Let us observe that q n verifies the following equation 

^q n = divu™ - J n (u n .Vq n ) - J n (q n dwu n ). 

According to the first part, (u n ) n ^n is u.b in so we can conclude that ^q n is 

u.b in L 2 {B~~ l ). Indeed we have: 

\\JJu n .V q n )\\ „ n , < ||u n || , jv \\q n \\ 

~JV_i JV JV 

And we recall that we use the fact that B 2 > 2 <^-> 2 . 

Let us prove now that Jrtf 1 is u.b in Ls{B~~^) + L^{B~~^-) and that <9 t r2 n is u.b in 

Lz{B^~z) (which gives the required result for u n in using the relation u n = — A _1 VcP — 

A^dM*"). 

Let us recall that: 

^d n = J n (u n ■ Vd n ) + J n k- l dw[A(v(p{l + q n )),U n ) + J n {K(<p(p(l + <f )))V(f )] 

+ vAd n + 6Aq n -R<f>* Aq n , 



-n n = J^curl [.4(^(1 + q n )), u n ) + J n (K(<p(p(l + g n ))V<f ))] + uAQ n . 
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Results of step one and an interpolation argument yield uniform bounds for u n in 

we infer in proceeding as for ^q n that: 

A n = J n (u n ■ VcP) + J n k- l d\v[A( V {p{l + q n )), u n ) + J n (K(<p(p(l + q n )))Vq n )] + *Ad" 

is u.b in Ls [B 2 2). 

Using the bounds for q n in L 2 {B^)r\L 00 {B^~ 1 ^), we get q n u.b in L 4 (B^~2~) in using 
proposition 12.61 We thus have J n (K(ip(p(l + q n ))Vq n u.b in L3(_Bt _ 2). 

N -1 , 4 , JVr 1 4iV3 

Using the bounds for u n in L°°(B~~ )nLs (iJT+2 ) we finally get A n is u.b in L3 (B~~2 ). 
To conclude </> * Ag™ is u.b in L i (B~~2) ) so jLd n is u.b in L~s (_B"2) + L 4 (B~~2). 

The case of goes along the same lines. As the terms corresponding to Aq n and 

* Acf 2 do not appear, we simply get c^Q n u.b in L3 (£? 2"~§). 

□ 

We can now turn to the proof of the existence of a solution and using Ascoli theorem to 
get strong convergence. We proceed similarly to the theorem of Aubin-Lions. 

Theorem 3.4 Let X a compact metric space and Y a complete metric space. Let A be 
an equicontinuous part of C(X,Y). Then we have the two equivalent proposition: 

1. A is relatively compact in C(X,Y) 

2. A(x) = {f(x); f £ A} is relatively compact in Y 

We need to localize because we have some result of compactness for the local Sobolev 
space. Let (Xp)peN be a sequence of Cq°(M. n ) cut-off functions supported in the ball 
B(0, p + 1) of R-^ and equal to 1 in a neighborhood of B(0,p). 

For any p S N, lemma [T] tells us that ((XpQ n , Xp un ))n<=N is uniformly equicontinuous 
in C(M + ; B~~ x x {B~~i) N ). In using Ascoli's theorem we just need to show that 
((XpQ n (t, •)> XpU 11 )^, -))nGN is relatively compact in B~~ l x (B~~2) N \/t £ [0,p]. 

~ JV -1 JV JV iV -1 

Let us observe now that the application u — > \ v u is compact from B 2 '2 = B 2 f]B 2 
into H~z~ l , and from B^^ 1 n B^~i into H~z~i. 

After we apply Ascoli's theorem to the family ((XpQ n i Xp nn ))neN on the time interval 
[0,p\. We then use Cantor's diagonal process. This finally provides us with a distribution 
(q,u) belonging to C(M + ; H^^ 1 x (H~~2) N ) and a subsequence (which we still denote 
by (q n , u n ) n £N such that, for all p £ N, we have: 

(x P q n ,X P u n ) -n^+oc (x P q,X P u) in C([0,p]; H^ 1 x (H^f) (3.34) 
This obviously entails that (q n ,u n ) tends to (q,u) in £>'(M+ x R N ). 

Coming back to the uniform estimates of step one, we moreover get that (q, u) belongs to: 

L^BT-Lf x (i? f+ 1 ) w )nL 00 ( J Bf^ 1 'f x (B^ +1 ) N ) 

and to C5(M+; x (C3(M + ; B ~2~i) N ). Obviously, we have the bounds provided of 

the firts step. 
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Let us now prove that (q,u) solves the system (SW), we first recall that (q n ,u n ) solves 
the following system: 

' d t q n + J n (u n ■ Vq n ) + divu n = - J n (q n divu n ) 
d t u n - i?Au n + 5Vq n - Rcj) * Vq n + J n {u n ■ Vu n ) + J n (K(ip(p(l + q n ))Vq n ) 

+ J„(^(l + 9 n )),/))=0 

The only problem is to pass to the limit in 2?'(R + x M. N ) in the non linear terms. This 
can be done by using the convergence results coming from the uniform estimates ()3.34p . 
As it is just a matter of doing tedious verifications, we show as a example the case of the 
term J n (K{<p(p(l + q n )))Vq n ) and J n (A(p(p{l + q n )),u n )). 
We decompose: 

J n (K(p(p(l + q n )))Vq n ) - K(p n )Vq n = J n (K(p(p(l + <f )))V<f ) - K(<p(p(l + q)))Vq. 

(Note that for n big enough, we have K(ip{p(\ + q n ))) = K(p n ) as we control ||/o n ||L°° 
and || — 1 1 2,°°). Next we have: 

J n (K(<p(p(l + q n )))Vq n ) - K(p(p(l + q)))Vq = J n A n + (J n - I)K{p{p{l + q)))Vq, 

where A n = K{p(p{l + q n )))Vq n - K(p{p(l + q)))Vq. 
We have then ( J n — I)K(p(p(l + q)))S7q tends to zero as n — ► +00 due to the property 

N 1 

of J n and the fact that K(p(p(l + q)))Vq belongs to L°°{B-- 1 ) L°°(L q ) for some 
q > 2. Choose V 6 C$°([0,T) x R N ) and p G C§°([0,T) x R N ) such that ip' = 1 on 
supp^, we have: 

I < (J n - I)K(p(p(l + q)))Vq, ip > I < ll^iY^l + g))^!!^^)!!^-/)^!!^, 

because L q loc <^-> L 2 oc and we conclude by the fact that ||(J n — I)il)\\ L 2 — ► as n tends to 
+00. 
Next: 



< J n A n ,i/j >=ll+j 



2 

n 1 M n 3 



with: 



4 =< (if (<?(p(l + <f ))) - i^(p(l + q))))Vq n , J n ^ >, 
J 2 =< K{p(p{l + <z)))V(<f - q), J n ^ > . 

We have then: 

^< 11^1^,11^(9"- dll^-ijM*-. 

Indeed we just use the fact that <p B^~ l and p H^~ l are embedded in L 2 . Next we 
conclude as we have seen that g n — ► ?l _ V 4. 00 q in Ci oc {H^ oc ). So we obtain: 

J*-„-+oo0 in D'((0,r)xl*). 

We proceed similarly for J 2 , indeed we have: 

In =< V(q n - q),p'div(K(p(p(l + q)))J n ip) > 
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and we have K(cp(p(l + q)))J n ip € L°°(B 2 ) so: 

II < \W\q n - g^^-ijUAWa + 9)))^V-ll Loo(B ^ } - 
We conclude then that: 

^^+co0 in v\(0,T*)xR N ). 

We concentrate us now on the term J n (A(<p{p(l + q n )),u n )). Let <p £ C§°(R+ x R N ) 
and p 6 N be such that suppc/? C [0,p] x l?(0,p). We use the decomposition for n big 
enough: 

tp J n A(<p(p(l + q n )),u n ) - yA{p, u) = tp' Xp A( V (p(l + q n )),X P {u n ~ «)) 

+ <p a(xp<p(pO- + 9 n )) - XpPO- + ?)),«). 

According to the uniform estimates and (13.34H . x p (u n — u) tends to in L 1 ([0,p]; H^ +l ) 
by interpolation so that the first term tends to in L (H 2 ) and we conclude for the 

- N 

second term in remarking that y- tends to ^ as p n in L°°(L°° n H~z). 
The other nonlinear terms can be treated in the same way. 

3.3 Proof of the uniqueness in the critical case 

Theorem 3.5 Let N > 2, and (qi,u{) and (q 2 ,u 2 ) t> e solutions of (SW) with the same 
data (qo,uo) on the time interval [0, T*). Assume that for i = 1,2: 

(q t ,u t ) £C([0,r),4,i) and u t £ {C([0,T*),B% 1 )nL} oc ([0,T*),B 2 NA )) N . 
There exists a constant a > depending only on N and physical constants such that if: 

Ikillfoo cri s < a, (3.35) 

then (qi,ux) = (q 2 , u 2 ) on [0, T*). 

N 

Let (qi,ui), (q 2 ,u 2 ) belong to with the same initial data, we set (5q,5u) = (q 2 — 
qi,u 2 — u\). We can then write the system (SW) as follows: 

d 

—5q + u 2 ■ V5q = Hi, 

dt (3.36) 
d 

—5u — vlS.5u = H 2 
< dt 

with: 

Hi = — div<5u — 5u ■ Vqi — 5qdivu 2 — qidivu, 

H 2 = —5V5q — R<fi* V5q — u 2 ■ V5u — 5u ■ Vui + A(q\, Su) + A(5q, u 2 ). 

Due to the term Su- Vq 1 in the right-hand side of the first equation, we loose one derivative 
when estimating 5q: one only gets bounds in L°°(B^ 1 ). 
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Now, the right hand-side of the second equation contains a term of type A(5q, U2) so that 
the loss of one derivative for 5q entails a loss of one derivative for 5u. Therefore, getting 
bounds in: 

C(R + ;B N \)nL 1 (R + ;B 1 NA ) 

for 5u is the best that one can hope. If enough regularity were available, we would not 
have to worry about this loss of derivative. But in the present case, the above heuristic 
fails because we have reached some limit cases for the product laws. Indeed, a term such 
as Su • V«i cannot be estimated properly: the product does not map x x 1 into 
Bft\ but in the somewhat larger space B^ 1 ^. At this point, we could try instead to get 
bounds for 5u in: 

c([o, r*) ; Brfj n LlM r*) ; b 1 ^), 

but we then have to face the lack of control on Su in L 1 (0, T;L°°) (because in con- 
trast with B\j 1 , the space Bj^ is not imbedded in L°°) so that we run into troubles 
when estimating 5u ■ Vq±. The key to that difficulty relies on the following logarithmic 
interpolation inequality (see the proposition 12.50 : 

+ l|U|l ^(^,oo) 



nil 7- 1 / R i \<||u||7i/ D i \ log e + 



and a well-known generalization of Gronwall the Osgood's lemma (see [2]) that we recall. 



Lemma 2 Let F be a measurable positive function and 7 a positive locally integrable 
function, each defined on the domain [to,t\]. Let \x : [0, +00) — > [0, +00) be a continuous 
nondecreasing function, with /z(0) = 0. Let a > 0, and assume that for all t € [toj^lL 



F(t) <a + [ ~f(s)fi(F(s))ds. 



If a > 0, then: 



ft r 1 ds 

-M(F(t))+M(a)< / j(s)ds, where Mix) = / — T . 

Jt Jx K s ) 

If a = and M(0) = +00, then F = 0. 
Proof of the theorem 13. 5t 

First step: in which space do we work? 

Let us observe first that in view of Sobolev embedding, £ C(M + ; L°°). Therefore, if a 
is small enough, by embedding and continuity we get: 

/ \ 1 1 
\qi{t,x)\ < - 

for x G R N and t in a small nontrivial time interval [0, T]. 

That observation will enable us to apply proposition [2?71 to the non-linear terms involving 
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ft- 

We shall further assume that T S (0, +00) has been chosen so small as to satisfy: 

C\\Vu 2 \\ L}r{Bki) <log2, (3.37) 

for some appropriate constant C whose meaning will be clear from the computations 
below. 

To begin with, we shall prove uniqueness on the time interval [0, T] by estimating (5q, 5u) 
in the following functional space: 

F T = L°°([Q,T];B° N>00 ) x (L"([0, T]; B^J n 

Indeed as explained below, in this space we can control the remainder because it is 
appropriate to the result of paraproduct. 

Why (5q,5u) is in Fy? 

Of course, we have to state that (5q, 5u) £ Ft, a fact which is not entirely obvious. We 
want now to show that (5q,5u) belongs to Ft- 

According to our assumption on (ftjUj), the estimates of paraproduct yield dtqi 6 
L T (B° Nl ). Therefore ft = ft — qo belongs to C2 ([0, T], B^ x ), which clearly entails by 
embedding Sq e C{[0,T),B° Noo ). 

Let tii = Ui — ul with ul solution to the following linear heat equation: 

d t u L - fJ.Au L = -SVq + kV(0 * q ), 
u L (0) = Uq. 

We obviously have (uj)o = and: 

d t Ui - vAiii = -Ui ■ Vui - <5Vft + RV((p * q,i) + A(pi, ft) + K(pi,Ui). 

The product and composition laws in Besov spaces insure that the right-hand side belongs 
to Lf^B^J (because B%^ x B^ ^ <-» BJ^^) thus to L^B^ 1 ^)) (for the last term, we 
use that ft G L^(B° Nl ). 
Now Proposition 12.91 implies that: 

Second step: Estimates on (5q, Su) 

Let us turn to estimate 5q. Proposition 12.81 combined with (|3.37p yields for t <T: 
)</ {\\ 5u -Vq\\ B o + \\6qdxvu 2 \\ B o +\\div5u\\ B o ))dr 

t V N,oo' N,oo iV,oo N,oo 

JV 

Estimate of type Bp j00 n L°° x i?p jOC) ^ -Bp,oo with s + ^ > enables us to get the 
following inequality: 

INL-Obo, )</ (INIbo ||div« 2 || B i nL"° + \\ Su \\b}, ni-( 1 + lki||Bi JdT, 

t ^ jV,oo' /(-, N,oo N,oo N,oo N,l / 
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whence, according to Gronwall inequality, to the embedding 1 Bjj ^ n L°° and to 
(|3.35|) we get: 

ll^ll^(^,oo) ~ ll 5u lli t l (Bj r>1 )( 1 + IkilL-^))- 

Making use of (|3.35p and proposition 12,41 we end up with: 

11 r 11 <u i, , ( , W%(*&„) + 

INLk^) < ||^||£x (Bkoo) iog^+ 



Remark that: 



with: 



||Hlzi (B o ) + ||HI £ i (B 2 ) <n*) = ^i(t) + v 2 (i) 

c ^ N .00' t \ N.oos 



Vi(t)= / (Ik(r)IU +||«i(r)|| B aJdr<+oo 



since L^(B° N1 ) ^ L\(B° N1 ) for finite t. 
We finally get: 



II^II^K^oo) ~ W^ ao ,l0g[e+ ynFT^ ), (3.38) 



with V non-decreasing bounded function of t G [0, +00). 

Let us now turn to the proof of estimates for Su. According to proposition 12.91 we 
have: 

ll^lk^oJ + II^H^^oo) ~ H U2 • V ^HZ t 1 (B^oo) + W A ( qi > Su) h(B*%) 

+ \\8u. ^MlI(b- n ^) + IW^IIzi^-ij + ll^Vftll^^, 

+ ||iC(g 2 )V^||r 1(B ! } . 

Let us assume that the a appearing in (|3.35p is small enough so that the second term in 

the right-hand side may be absorbed by the left-hand side. ||u2||r2/ R i n tends to when 

t y n,i' 

t goes to 0, so if we choose T small enough, the first term may also be absorbed. Using 
interpolation of proposition 12.41 we obtain for all t £ [0, T], 

H^lk^oo) + II^H^(^,oo) ~ I DMIb^JMIb^ + (1 + IMI^INIb^J^ 

Let us now plug (|3.38p in the above inequality. Denoting: 

X(t) = \\Su\\ L?{B -^ ) + \\Su\\ Ll ^ ) 

we get for t < T, 

X(t) < £ (1 + V'(r))X{r) \og(e + ^)dr 
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As: 

dr 

V(T) I CXZD , 

r logfe H — 



V , €L 1 (0,T) and f 
Jo 



Osgood's lemma (see lemma [2|) implies that X = on [0, T], whence also Sq = 0. 
Standard arguments of connexity then yield uniqueness on the whole interval [0, +oo). 



4 Proof of theorem 11.2 

We will proceed similarly to the proof of the theorem 11.11 To begin with, let us observe 
that under the definition 11.21 system reads: 

' dtq + u ■ Vu = H 

< d t u - uAu = -K (4.39) 
k (q,u) t =o = {qo,u ) 

with: 

H = -(l +q)divu, 

K = G — u- Vu + (P' (p) + n)Vp - Kcp * Vp. 

As previously we can build approximate smooth solutions (q n ,u n ) of (|4.40p in studying 
the Korteweg system with a capillarity coefficient K n = —. It is convenient to split (q n , u n ) 
into the solution of a linear system with initial data (q$, Uq), and the discrepancy to that 
solution. More precisely we denote by (g£, u^) the solution of the linearized pressure-less 
system on the intervall [0, T]: 

' d t q n L + div< = 0, 

< d t u n L - pZAul — (A + fl)Vdivul = 0, 

>2,<)a=o = (M), 

with: 

(% n ^o) = (E A ^o ; E A < n °)- 

\l\<n \l\<n 

We set: 

(q n ,u n ) = (q n L + q n ,< + u n ) 
We can state now that (q n ,u n ) verifies the following linear system: 

d t q n + divn n = F n , 

d t u n - flAu n - (A + /2)Vdivn n - -VAq n = G n , (4.40) 

n 

(<r,u n ) /t=0 = (o,o) 

with: 

F n = -q n divu n , 

G n = -u n ■ Vu n + A{ P n , u n ) - K(p n )Vq n + (P' {p) + K )Vp n - K(f) * Vp n . 

We want show that such solution (q n ,u n ) exists, in this goal we recall some theorem by 
R. Danchin and B. Desjardins in [19j . 
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Theorem 4.6 Let p G [l,+oo[. Then there exists r\ > such that if qo G -B ? , uq G 



AT 



(Bl^ 1 )^ and: 



ko w < 77, 



7-V- 



i/ien i/iere exists T > suc/t £/tai system /ji4-40\ ) has a unique solution (q, u) in E 

In fact, we can extend this result to the case where the viscosity coefficients are variable. 
And we can show that the uniform estimates are independent of the capillarity coefficient. 
We can obtain the following result on our solution [q n ,u n ) n ^. 

N_ 

Theorem 4.7 Let p G [l,+oo[. Then there exists r\ > such that if qft G B p , Uq G 
(BT ) N and: 

H\\ B f<v, 



then there exists T > such that system ^4-4$ has a unique solution (q n ,u n ) in E^ n 
and (q n ,u n ) are uniformly bounded in Fj,. 

Proof: 

We recall a proposition (see [19J) on the following linearized pressure-less system: 

dtq + divn = F, 

dtu — pAu — (A + p)Vdivu — nVAq = G. 



Proposition 4.12 Let s G M, p G [1, +00], 1 < pi < +00 and T g]0, +00]. // (go, ito) G 
B^ x (B^ 1 ) 1 * and (F,G) G L^^Bp" 2 ^ x (Bp~ 3+1 *) N ), then the above linear system 

has a unique solution (q,u) G C T (B° x (Bp" 1 ) 1 *) n L^(Bp +JI x (B v 1+Tl ) N ). 
Moreover for all p G [pi,+oo], there exists a constant C depending only on p, X, R, p, 
pi and N such that the following inequality holds: 

\\q\\~ p s +| + \\ u \\~ p s -i+f) ^ c(lkol|-B£ + ||«o|| B ;-i 

Lrjp {Bp ) Lrjp (Bp 



+ \\F\\_ K _^2_ +IIGI 



Remark 2 More precisely we have: 

Ml¥(B») + K ll9llzi.(B^+ 2 ) ^ c (W<1o\\b° + ||«o|Ibj-i + ll-^llz^^i) + HGIIl^- 1 ))- 
Uniform Estimates for (q n ,u n ) ne n 

The existence of solutions in the case of general viscosity coefficient follow the same line 
as the proof in [19]. It suffices to solve a problem of fixed point. 
Denoting by V(t) the semi-group generated by system (|4.4ip . we have: 

{q%ul){t) = V(tM,u%). 
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Let us define: 

<P ql , ul (q n , u n ) = f V(t - s)(F(ql + q n ,< + u n )(s), G(q n L + q n ,< + u n )(s)) ds. 
Jo 

where we have set: 

F{q,u) = -div(gu), 

G(q, u) = -u-Vu + +A{p, u) - K(p)Vq + (P (p) + n)Vp - n<p * Vp. 

In order to prove the existence part of the theorem l4.71 ther's just have to show that </>g£,u£ 
has a fixed point in Fj,. Since Fj, is a Banach space, we can prove that <Aj£,u£ satisfies the 
hypothesis of Picard's theorem in a ball (B(0, R) of F^ for sufficiently small R. Moreover 
R depends only of mathematical constants. We can find a time T independent of n such 
that for all initial data verifying H^qH n < -f , we have existence of solution (q n ,u n at 

less on the interval (0, T). The end of proof consists to verify that (q n , u n ) is uniform in 
Fj,, it suffices to use the proposition 14.121 

We now have builded approximated solution (g n , u n ) of system (SW) and we can conclude 
in using technics of compactness. 



2) Existence of a solution 

The existence of a solution stems from compactness properties for the sequence (q n , u n ) n ^fq 
an we want use some result of type Ascoli as in the proof of theorem 11.11 

Lemma 3 The sequence (dtif 1 ,dtu n ) n ^ is uniformly bounded in: 
L 2 (0,T;£p'~ _1 ) x (L^T;^-- 1 '-- 2 ))^, 

for some a > 1. 
Proof: 

Throughout the proof, we will extensively use that Lj(Bp) <^-> Lj(Bp). The notation u.b 
will stand for uniformly bounded. 
We have: 

t q n = -u n ■ Vq n - (1 + q n )divu n , 

1 (4-42) 
d t u n = -u n ■ Vu n - q n A(p n ,u n ) - K(q n )Vq n + - VA«f. 

n 

~K iY._i 

We start with show that dtqf 1 is u.b in L 2 (0, T; B v v ' v ). 

9 — — -1 

Since u n is u.b in L^(Bp) and X7q n is u.b in L^?(i?p ), then u n ■ Vq n is u.b in 

~jy_ JV_ 1 

Lj(B p ' p ). Similar arguments enable us to conclude for the term (1 + q n )divu n which 

^JV N_ — \ K ——1 

is u.b in h\,{B v v ' v ) because q n is u.b in ^(Bp ) and dfvu n is u.b in h\(B v v ). 



Let us now study dtu n+ . According to step one and to the definition of u^, the term 

jv_ 2 — — 1 — — 1 

_ _ . _ - ,v 

'P 



Au n+l is u.b in L 2 (B P P ). Since u n is u.b in L°°{B P P ) and Vn n is u.b in L 2 (B P P 
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——2 ~2L — i K—2 

so u n ■ Vu n is u.b in L 2 (B P P ) thus in L 2 (B P P ' p ). 

JV 

Moreover we have q n is u.b in L°°(Bp ) and g n is u.b in L°°, so by proposition 12.71 

N_i „2£_l JV_2 

VK (q n ) is u.b in L OD (B p p ) thus in L 2 (B P P ' p ). This concludes the lemma. □ 

Now, let us turn to the proof of the existence of a solution for the system (SW). We want 
now use some results of type Ascoli to conclude in use the properties of compactness of 
the lemma El 

According lemmaEJ (q n ,u n ) n ^ is u.b in: 

_jv JV_ 1 . ^jX_i iX_2 

CmO,T];B p p '» )x(C 1 --([0,T];B p " p )f , 

thus is uniformly equicontinuous in C(([0, T]; Bp ' p ) x {Bp ' p ) N ). On the other 
hand we have the following result of compactness, for any <j) E C^(R N ), s £ M., 5 > the 
application u — > <pu is compact from B^ to Bp S ~ . Applying Ascoli's theorem, we infer 
that up to an extraction (q n ,u n ) n( zfq converges in 2?'([0,T] x M N ) to a limit (q,u) which 
belongs to: 

__jy_ n__i , ,^JY._i K—2 

CH[0,T];B p p,p )x(C 1 ~([0,T];B p * ' p )) N 

Let (q, u) = (q,u) + (qo,UL). Using again uniform estimates of step one and proceeding 
as, we gather that (q, u) solves (SW) and belongs to: 

p + L?(B p p,p )x(L 1 T (B p p )PiL?(B p p )) N . 

Applying proposition, we get the continuity results: 

__jy_ JV_-l JV j 

p-peC([0,T],Bp p ' p ), ueC([0,T],B p p ). 

5 Proof of theorem 11.31 

In this section, we consider the case when the initial density belongs to p + B p ' ? and 
satisfies < p < po- We consider again the study of a approximate sequence verifying 
(|4.40p and we proceed as previously. 

Construction of approximate solutions 

We consider the following system: 

d t q n+1 + u n ■ Vq n+1 = q n dwu n , 
8 t u n+l + u n ■ Vu n+l - A(u n+1 ) + Vq n+1 = G n , 

with: 

G n = A(u n ) - A(p n , u n ) + K(p n )Vq n . 



(5.43) 
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Uniform estimates for (g ra , u™)„ e N 

Let us show that the sequence (q n ,u n ) n( z^ is uniformly bounded in provided that T and 
77 have been chosen small enough. 

Let us remark first that according to proposition, there exists a universal condition K 
such that for all n € N, we have: 

K||~ jv <U , 



with C/q = ^ll n o|| ^ , 



<n zl r^^^s: 29( ^ 1)||A ^ 01 



Suppose that II go II jv < n for a small constant 77 and let Co = 1 + II go 1 1 jv • According to 

Bp Bp 

, we can choose a positive time T such that the following property holds for all n G N: 

2 i-i 

IbrlL jv < Uq, and llnrlL n,, <V r U n r 
Let us now show by induction that the following estimates are satisfied: 

||g n ||~ K <^q 

ll« n H~ N , + \\u n \\„ N 1-1 < 77 
L™{BT 1 ) U L 1 T (BT + 1 ) ~ 

From now, we suppose that 7/ < (2C±) where C\ is the norm of the injection B p L°°. 
This ensures us that the following inequality is satisfied: 

According to proposition, we have: 

C||«"ll , n +1 

\\q n+1 \\~ Jv <exp l t^ bp )(||g || n + \\F n \\ m) 
L™(BT) ~ " yu ^' u BT "L^,(BT) ; 

Moreover we have by Proposition: 

C*|K|| N +1 

||u n ||_ jv <exp l t^ bv )(|| go || at +11^11 n) 

L™(B~) V B~ L\,(BTy 

whence according to proposition: 

\\F n \\ , jv <T 1_ n(l + ||g n || jv )||div« n || n jv < (1 + 170)5*7, 

"4(BP) "l°?(Bp) L r ^(BP) 



and we get finally: 



|g n |U jv <{l + U Q ) 1 2r 1 e^ c{l+u ^ ri 

L™(B~) ~ 
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Obviously if i] has been chosen small enough then q n+1 satisfies the estimate in (J-'n+i)- 
Applying proposition to the second equation of yields: 

||n n+1 |U n +\\u n+1 \\ N+1 <\\u n -Vu n \\ n , + \\K {q n )Vq n \\ 



+ 'I" 4 *"" (B *-v 



We now use proposition 12.71 as in the proof of ?? to conclude. 



Existence of a solution 

We can now easily show that (q n ,u n ) is a Cauchy sequel in our space Ft of uniqueness 
and so (q n ,u n ) -> (q,u) in F T . 

It rests to verify by compactness that (q,u) is a solution of the system (SW). □ 
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